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hyperkagome lattice: Comparison with Na4lr308. 
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We develop high temperature series expansions for In^ and the uniform structure factor of the 
spin-half Heisenberg model on the hyperkagome lattice to order (3^^ . These expansions are used to 
calculate the uniform susceptibility (x), the entropy (S), and the heat capacity (C) of the model 
as a function of temperature. Series extrapolations of the expansions converge well down to a 
temperature of approximately J/4. A comparison with the experimental data for Na4lr308 shows 
that its magnetic susceptibility is reasonably well described by the model with an exchange constant 
J fti 300iC, but there are also additional smaller terms present in the system. The specific heat of 
the model has two peaks. The lower temperature peak, which is just below our range of convergence 
contains about 40 percent of the total entropy. Despite being a 3-dimensional lattice, this model 
shares many features with the kagome lattice Heisenberg model and the material must be considered 
a strong candidate for a quantum spin-liquid. 



PACS numbers: 74.70.-b,75.10. Jm,75.40.Gb,75.30.Ds 



INTRODUCTION 

Despite many experimental and theoretical studies, the 
spin-half kagome lattice Heisenberg Model remains in- 
completely understood and a subject of intense ongoing 
research. [l|-[3| The kagome lattice is a two-dimensional 
lattice of corner sharing triangles, where classical anti- 
ferromagnetic models lead to exponential ground state 
degeneracy. Whether quantum fluctuations lift this de- 
generacy and lead to an ordered state, or select a valence- 
bond order, or lead to a resonating valence-bond quan- 
tum spin-liquid possibly with topological order and frac- 
tional excitations remains a topic of great interest. Nu- 
merical evidence shows that many distinct singlet states 
are very close to each other in energy, (i- 13 1 

The hyperkagome lattice is a three-dimensional lattice 
of corner sharing triangles, where similar ground state 
degeneracy at the classical level and possible compet- 
ing phases with nearly the same energy in presence of 
strong quantum fluctuations remains valid. Experimen- 
tal discovery 14] of the material Na4lr308, where the Ir 
ions form a hyperkagome lattice, with strong signatures 
of a gapless quantum spin-liquid phase has led to much 
excitement. Theoretical proposals for this system also 
include various quantum spin-liquids as well as valence 
bond order, isl- isj 

In this paper, we develop high temperature series ex- 
pansions for In Z and the uniform structure factor of the 
spin-half hyperkagome lattice Heisenberg model to order 
/3^^. These series expansions are used to calculate the 
uniform susceptibility, the entropy and the heat capac- 
ity of the model. Pade and d-log Pade approximants 
show that series extrapolations converge well down to 



a temperature of approximately J/4. At this tempera- 
ture 40 percent of the total entropy of the system is still 
present. The heat capacity, which has a short-range or- 
der peak around T = 2J/3 shows hints of another rise 
below T = J/4, leading to a sharp peak in C/T. These 
results are reminiscent of similar results on the kagome 
lattice. E3421 



Although the susceptibility data can be well fit by a 
Curie- Weiss law over a wide temperature range, the effec- 
tive Curie- Weiss constant keeps changing with the region 
of fit. [is! Comparison with the experimental susceptibil- 
ity data shows that the material is described by a J value 



of approximately 300i<r.[15|,[16| This J value is consistent 



with more recent high-field measurements. [2^ There are 
also clear deviations below T = J/2, where the series ex- 
trapolation are still well convergent. These presumably 
arise from smaller exchange terms induced by spin-orbit 
couplings, such as the Dzyaloshinski-Moria interactions, 
which are not included in our theoretical studies. 2^ 26 1 
Similar anisotropy terms are also known to be present in 
the kagome antiferromagnets.[3, 22, 27| 

Experimental determination of the magnetic heat ca- 
pacity requires subtraction of non-magnetic contribu- 
tions, which becomes increasingly uncertain with increas- 
ing temperature. Our study can not shed light on the 
low temperature power-laws observed in the material. At 
temperatures above J/3 which translates to ~ lOOK for 
the material, where our results are well converged, there 
is a clear discrepancy between the experiments and the 
model calculations. We show that these discrepancies 
can be reconciled if one assumes some additional har- 
monic modes below 300K that have not been fully elim- 
inated by the subtraction procedure in the experimental 
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paper. 14 1 



MODEL AND SERIES EXPANSIONS 

We study the nearest-neighbor spin-half Heisenberg 
model 



7/ — J ^ ^ Si • Sj 



(1) 



Here, the sum runs over each nearest-neighbor pair of 
spins on the hyperkagome lattice. We develop high tem- 
perature series expansions for hi Z and the uniform struc- 
ture factor, which is also Tx, where x is the uniform 
susceptibility. These quantities are given by 



In Z = In Tr exp — 131-1, 



and 



TX = Y.< ^^^^ 



> 



(2) 



(3) 



In 2 

O.OOOOOOOOOOOOD+00 
0.300000000000D-I-01 
O.OOOOOOOOOOOOD-hOO 
-0.850000000000D-f01 
O.OOOOOOOOOOOOD-fOO 
0.465333333333D4-02 
0.213333333333D-f01 
-0.312973809524D-f03 
-0.354539682540D-f02 
0.235000994709D-h04 
0.467922116402D+03 
-0.188903919384D-h05 
-0.570519634440D+04 
0.158860661564D-H06 
0.668440957641D-h05 
-0.137796110009D-f07 



0.25 

0.100000000000D4-01 
0.200000000000D+01 
O.lOOOOOOOOOOOD+01 
-0.666666666667D+00 
0.134666666667D-f02 
0.352888888889D+02 
-0.699206349206D+02 
-0.237619047619D-I-03 
0.745729100529D-H03 
0.236886624339D-H04 
-0.685090572391D+04 
-0.231225320571D+05 
0.656040985160D4-05 
0.236506418283D-f06 
-0.623070199393D-f06 
-0.244735886559D+07 



TABLE I: Series expansion coefficients for InZ and uni- 
form structure factor for the hyperkagome lattice Heisenberg 
model. 



Here the sum runs over all the spins of the lattice, and 
the angular brackets denote a thermal expectation value. 
Details of the method for calculating the series expan- 
sions can be found in the literature. 1 28 Let us write 



InZ 



where N is number of sites, and. 



(4) 



(5) 



The coefficients a„ and 6„ upto n = 16 are given in Table 
1. Note that the series coefficients are valid for both 
the ferromagnetic and the antiferromagnetic couplings. 
The factors of (—1)" are needed for the antiferromagnetic 
model, where J > 0. 



SERIES ANALYSIS AND RESULTS 




log(T/J) 



FIG. 1: Variation of the efi'ective Curie- Weiss constant with 
temperature for the Heisenberg Antiferromagnet on the hy- 
perkagome lattice. 



We begin with an analysis of the uniform susceptibil- 
ity series. To convert the susceptibility into emu/mole, 
we need to multiply Tx by a factor of Naq"^ fJ-%/4:kB — 
0.0938(7^. At high temperatures the susceptibility can be 
fit to a Curie- Weiss law, whose asymptotic value for our 
model is — J. However, as discussed by Zheng et 
al, 23] the effective Curie- Weiss parameter obtained by 
fitting the inverse susceptibility versus temperature data 
to a straight line slowly changes with temperature. It 
can be defined around a temperature T as 



A plot of this effective Curie- Weiss constant as a func- 
tion of log(T/J) is shown in Fig. 1. We see that even at 
temperatures of lOJ the effective Curie- Weiss constant 
is 10 percent off its asymptotic value. At temperatures 
around and below J the effective Curie- Weiss constant 
is larger than 2 J. This means that the exchange con- 
stant must be a factor of 2 or more smaller than de- 
termined by Okamoto et al.[l3l Indeed, J = 300K and 
g = 1.99 leads to an excellent fit of the experimental data 
above T = 150^. Since the g-factor is not independently 
known, a finer fit is not useful. We note that the series 
extrapolation converges well down to about J/4 and be- 
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FIG. 2: A comparison of the uniform susceptibility data of 
Na4lr308 versus temperature with the Heisenberg model for 
J = 300K and g = 1.99. 



FIG. 3: The molar heat capacity for the hyperkagome lattice 
Heisenberg model with J = 300K. The [N+2/N] pade ap- 
proximant ensures a behavior at low temperatures. But, 
it has a 40 percent missing entropy. The d-log Pade approxi- 
mants imply a second peak below T = lOOK. 



low J/2, there are clear deviations between the exper- 
imental data and the Heisenberg model. We attribute 
these to spin-orbit couplings, which must be present in 
the material. "265 Nevertheless, the susceptibility compar- 
ison shows that the nearest-neighbor Heisenberg model 
is a good starting model for the material, [iij] All other 
terms are presumably an order of magnitude smaller. 

For the rest of the study, we will fix the J value to 
SOOi^T. Series extrapolations for the heat capacity are 
shown in Fig. 3. There is a peak associated with short- 
range order around a temperature of 200K. This peak is 
well captured in our calculation. However, the extrapola- 
tions do not converge well below the peak. Experiments 
show a heat capacity at low temperatures. If we use 
an [N-I-2/N] pade to represent the heat capcity, it en- 
sures a dependence. All such extrapolations are close 
to the [6/4] approximant shown. However, they can not 
be right. The integration of the heat capacity implies 
a large missing entropy of order 40 percent. A similar 
problem was observed in extrapolating the high temper- 
ature series for the kagome lattice Heisenberg model. 
If instead we use d-log Fade approximants to extrapolate 
the series, we find that the majority of the approximants 
show a second rise in the heat capacity at a tempera- 
ture below lOOK. However, once again the convergence 
breaks down as the second peak arises. The method de- 
veloped by Misguich and Bcrnu[21] may be helpful in 
this regard. However, their method is based on taking 
advantage of our knowledge of the ground state and the 
low temperature properties in the extrapolation of the 
high temperature expansions. Since we know even less 
about the ground state properties of the hyperkagome 



lattice than we do about the kagome lattice, it is unclear 
how useful it would be at this stage. 

A comparison of the experimental C/T data with the 
theory is shown in Fig. 4. We clearly cannot address 
the low temperature power-laws, which occur at temper- 
atures well below our range of convergence. The discrep- 
ancy at higher temperatures can possibly be attributed 
to an incomplete subtraction. Experimental determina- 
tion of the magnetic specific heat requires subtraction of 
phonon and other contributions. When the J values are 
comparable to room temperature, this becomes difficult 
as the phonon contributions are not simply power-laws. 
For this material, the non-magnetic Na4Sn308 was used 
to subtract the phonons. But such a subtraction need not 
be accurate above some temperature. In fact, as shown in 
Fig. 5, the entropy obtained from the experimental data 
exceeds the theoretical curves above 70 K and exceeds 
In 2 before 200 K. In our calculations, the entropy is well 
converged down to about 70 K. It is also much less sensi- 
tive upto these temperatures to extrapolation methods, 
and to addition of other smaller exchange constants. 

We believe, there can be two reasons for this discrep- 
ancy. First, the subtraction is incomplete. Indeed, if we 
further subtract harmonic oscillator modes with a den- 
sity of states of about 0.25 moles at a single frequency of 
200K that can bring the data into good agreement with 
experiments, as shown in Fig. 4. A second possibility 
is that the materials have additional entropy associated 
with itinerant degrees of freedom because they are not 
too far from a metal insulator transition. This issue can 
presumably be addressed by further experiments. 
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the heat capacity and the entropy down to a tempera- 
ture of J/4. Comparison with the experimental data for 
Na4lr308 suggests that the Heisenberg model is a good 
starting model for the material with an exchange con- 
stant J w 300K. However, we find that there must be 
additional terms also present in the material. The ob- 
served low temperature power-laws can not be addressed 
by our calculations. We have suggested possible resolu- 
tions to discrepancies in higher temperature heat capac- 
ity and entropy data for the material. 

We would like to thank Dr. Y. Okamoto and Dr. H. 
Takagi for many fruitful communications and for send- 
ing us their experimental data. We are grateful for the 
computing resources provided by the Australian Partner- 
ship for Advanced Computing (APAC) National Facility. 
This work is supported in part by NSF grant number 
DMR-1004231. 



FIG. 4: Comparison of C/T between theory and experiments. 
The line with filled circles is the experimental data from 
Refill. The solid lines show representative behavior from se- 
ries extrapolations. The dashed line is experimental data from 
which additional subtractions corresponding to a heat capac- 
ity of 0.25 moles of phonons at a single frequency {200K) has 
been carried out. 
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FIG. 5: The entropy function calculated for the hyperkagome 
lattice Heisenberg model and the experimental material. 



CONCLUSIONS 

In this paper, we have developed high temperature 
series expansions for various thermodynamic properties 
of the Heisenberg model on the hyperkagome lattice. 
These expansions enable us to compute various thermo- 
dynamic properties, such as the uniform susceptibility. 
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